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Vorrichtung zum Erzeugen tomografischer Bifder 

@ Vorrichtung zum Erzeugen tomografischer Bilder, die die 
Verteilung der elektrischen Impedanz einer Korperschicht in 
ejner Ebene (x, y) eines Untersuchungsobjektes (17) darstei- 
ten, durch Messung von elektrischen Spannungen an der 
Oberflache des Untersuchungsobiektes (17) bai unters^ied- 
lichen Prequenzan, mit: 

— Elektroden (1 — 16). die am Umfang des Untersuchungs- 
objektes (17) enttang des Randes der Korperschicht anieg- 
bar sind, 

— einer Ansteuer- und MeQeinheit (13). die Wechselstrom 
von einem Stromgenarator (20) wenigstens zwei der Elektro- 
den (t, 2] zufOhrt, ferner zwischen jeweils zwef der restlfchen 
Elektroden infotge des Wechselstroms vom Stromgenarator 
auftretenda Spannungsdifferenzen mi8t und entsprechende 
MaSsignala lief art, 

— einer Elektroden umschalteinheit, die fortlaufend andere 
Elektroden mit dam Stromgenarator (20) verbindet, und 

— einer Auswerteeinheit (22) zum Verarbeiten der Me&al-. 
gnale fOr e!ne Anzeige des tomografischen BHdes auf einar 
Anzsigeeinhett (23), wobei die Auswerteeinheit (22) iewells 
eine den jewetligen gemessenen Potentialdifferenzen zuga* 
ordnete komplexe Impedanz oder Admittanz mit Real- und 
JmaginSrteii oder Betrag und Phase ermftteit, 
dadurch gekennzeichnat, 

da& die Frequenzen so gewihit sind, dad sie wenigstens 
etn Stuck des Dispersionsgebiets des Untersuchungsobjekts 
einschlie&en, und 

da& die Auswerteeinheit (22) den Imaginarteil oder statt 
der Phase das Verhattnis von Imaginartei) zu Raatteii durch 
die Froqaenz dividiert, so dafi auf der Anzeigeeinhelt (23) ein 



solches tomografisches Bitd angezeigt wird, das aus Bild- 
punkten aufgebaut tst, die durch eine Ruckprojektion auf 
Grund von zwei Messungen bei den gewahlten Frequenzen 
aus mindestens den Imaginarteil und/oder die Phase enthal- 
tanden Impedanz- oder Admittanzwerten gewonnen sind. 
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' Beschreibung * 

Die Erfindung geht aus von einer Vorrichtung zura Erzeugen tomograf ischer Bilder nach dem Oberbegriff des 
Patentanspruches 1. 

5 Die elektrische Impedanz von bioiogischem Gewebe ist schon seit langem Gegenstand lunfangreicher Unter- 
suchungen (vgl. z. B, Schwan, RP.: Electrical Properties of Tissue and Cell Suspensions. In: Advances in Biologi- 
cal and Medical Physics, J.H. Lawrence and CA. Tobias (eds.), VoL 5, New York (1957) 148ff, and Pethig, R.: 
Dielectric and Electronic Properties of Biological Materials. Wiley, Chichester (1 979)). 
Auch in der Medizin hat die Auswertung der Impedanz schon seit Jahrzehnten beispielsweise bei nichtinvasi- 

10 ven Methoden zur Bestimmung des Herz-Zeit-Volumens Anwendung gefunden (vgL z. B. Kubicek, W.G.: Deve- 
lopment and Evaluation of an Impedance Cardiac Output System. Aerosp. Med. 37 (1966) 1208 — 1212). 

Da verschiedenartige Gewebe sich in ihrer Impedanz unterscheiden, wird diese Eigenschaft seit einigen 
Jahren auch zur Erzeugung von Bildem benutzt, d. h. die Verteilung der Impedanz in einer Korperschicht n^rd 
zur Abbildung gebracht (vgL z. B. Barber, DSl, BJL Brown: Recent Developments in Appfied Potential Tomo- 

15 graphy-APT* Information Processing In Medical Imaging ed S. Bacherach^ Martinus Ntjhoff, Dordrecht (1986) 
106—121 und Webster, J.Gj Electrical Impedance Tomography, Verlag Adam Hilger, Bristol und New York 
1990), hierzu werden iiber <Ke Oberflache eiiies Untersuchungsobjektes diesem Wechselstrome eingeprSigt und 
an dessen Oberflache Spannungen gemessen, die durch die Strome und durch die Widerstandsverteflung im 
Kdrperinneren des Untersuchungsobjektes bedingt sind (vgl z. B, GB-A-2 119 520, GB-A-2 160323 und US- 

20 A-4 617 939), 

Eine prinzipieUe Schwierigkeh bei der Verwendung des elektrischen Stromes als Mittel zur Abbildimg liegt im 
Gegensatz zu Rontgenstrahlen darin, daB der Weg des Stromes durch das Korperinnere im allgemeinen nicht 
bekannt ist Die Stromdichteverteilung hangt namlich bei gegebener Stromeinspeisung von der ortlichen Verte?i- 
lung der komplexen Leitfahigkeit ab, die aber gerade in Erfahrung gebracht werden sollte. Trotzdem wurde 

25 nicht ohne Erfolg versucht, mit vereinfachenden Annahmen geeignete Aigorithmen zu entwickein, die, einge- 
setzt in fur diesen Zweck konstruierte MeBgerate, Bilder der Widerstandsverteilung im fCqrperinneren erzeugen 
kdnnen (vgL obige GB- und US-Schrif tenX 

Beim derzeitigen Stand der Technik werden Bilder aus dem Realteil oder dem Betrag der Impedanz erzeugt 
und zwar uberwiegend bei einer festen Frequenz des eingespeisten Stromes in der GrdBenordiiiing. von 50 klfa, 

30 fur die technische Grunde sprechen. Die imaginare Komponente (oder auch die Phase) wurde bisher nicht zur 
Abbildung herangezogen, weil sie im VerhSltnis zura Realteil klein erscheint und wenig Information enthalten 
soli (vgL Brown, B.H., AJD. Seagar: The Sheffield Data Collection System. Clin. Phys. PhysioL Meas. 8 SuppL 
A(1987)9l--97. 

In IEEE Transactions on Biomedical Engeneering, Band 40, Nr. 1, Januar 1993, 5. 29—34 wird ein elektrisches 
35 Impedanztomographieverfahren vorgestelit, bei dem die Impedanzverteilung eines Objekts durch Einpragen 
eines Wechselstroms und Messen einer korrespondierenden Spannung ermittelt wird Die hierbei verwendete 
Vorrichtung weist die ira Oberbegriff des Anspnichs 1 angegebenen Merkmale auf. Die MeBelektroden werden 
kreisformig um das Objekt angeordnet, und ihr MeBsignal wird bei einer Frequenz von 50 kHz ausgewertet 
Dabei werden die Realteile der an den Elektroden anliegenden Spannungen gemessen, und somit lediglich die 
40 ohmschen Anteile der zu messenden Impedanzen verarbeitet Nach dein Finite-Elemente-Verf ahren zusammen 
mit dem Newton-Raphson-Algorithmus werden den Impedanzmefiwerten entsprechende Bildpunkte ziigeord- 
net. 

Die Druckschrift WO 92/19153 zeigt eine Vorrichtung zur Bestimmung der komplexen Impedanz und der 
entsprechenden physikalischen Charakteristik eines MeBobjekts in Abhangigkeit von den verwendeten Pre- 
45 quenzen. Mit Hilfe von zwei Elektroden wird der eingespeiste Strom und die korrespondierende Spannung am 
Objekt gemessen. Die nach einer Kreuzkorrelation und Fahung resultierende komplexe Impedanz wird fiir die 
jeweiiige MeOstrecke an einer Sffemanzetge dargestellt 

Aus J. G. Webster: "Electrical Impedance Tomography Verlag Adam Hilger, 1990, Seiten 8 bis 14 und 1 14 bis 
117, sind Vorrichtungen zum Erzeugen tomographischer Bilder der eingangs genannten Art bekannt Zwar ist 
50 dabei auch angegeben, daB die Dielektrizitatskonstante sich frequenzabhangig andert Jedoch wild nur der 
Realteil der Impedanz zur Auswertung herangezogen* 

Es ist Aufgabe der vorliegenden Erfindung, eine Vorrichtung zur Erzeugung tomografischer Bilder zu schaf- 
fen, die sich durch hohe Qualitat und groBe Aufldsung auszeichnea 

Diese Aufgabe wird rait dem Gegenstand des Patentanspnichs 1 gel5st. 
55 Vorteilhafte Ausfiihrungsarten der Erfindung sind in den Unteranspruchen angegebea 

Die beanspruchte Vorrichtung, die auBer dem Realteil bzw. Betrag der Impedanz auch den Imaginarteil bzw» 
die Phase zur Abbildung auswertet und auBerdem den FFequehzbereich von 50 kHz bis hinab zu 10 Hz ausdehnt, 
erlaubt die Erzeugung Zustands- oder Frequenz-differenter Bilder. Es werden beispielsweise 16 Oberflachen- 
Elektroden verwendet, die zur Stromeinspeisung oder Spannungsmessung angesteuert sind, um die Mdgiichkeit 
60 zu haben, verschiedene Aigorithmen zur Bildverarbeitung einsetzen zu kdnnen; bevorzugt wird aber ein Algo- 
rithmus der "Backprojection". 

Die theoretischen Grundlagen fur die Impedanztomografie folgen aus der elektromagnetischen Feldtheorie. 
Das Ohmsche Gesetz und die Maxwell-Gleichungen beschreiben bekanntiich fur ebene Felder den Zusammen- 
hang zwischen elektrischem Feld E, magnetischem Feld H, elektrischem FluS D, magnetischera FluB B, einge- 
65 pragter Stromdichte J und einer Leitfahigkeit IC Die quasistationare Naherung wird angenommen, um das 
elektrische Feld als Gradient eines Potentials darzustellen. 

Die Maxwellgleichungen fCir isotrope, Ibeare Leiter lauten nach dem Obergang zur quasistationilren NShe- 
rung (ft, 6, 3 = komplexe Amplituden) 
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S-D = 0, axH = ^icE+j) l.Maxwellgleichimg (1) 



c 



S-B = 0, SxE = 0 Induktionsgesetz (2) 

mit der komplexen Leitfahigkeit 



10 



^"^^^4ii' a = Ohiusche Leitfahigkeit, (3) 

6) = Kreisfrequenz 

15 

Aus dem Induktionsgesetz (2) kann die komplexe Amplitude der elektrischen Feidstirke £ durch das skaiare 
Potential <Dgema£ 

(4) 

lokal integriert werdea Zusammen mit der Kontinuitatsgleichung 

a-(KE) = -d-j (5) 

erbalt man somlt fur das Potential <P die lineare Diff erentialgleichting 

30 

a.(Tt3o) = K3^*+(a-K)(8«>) = a-j (6) 

Wird Gleichung (6) durch h dividiert und mit xo als homogener Leitfahigkeitsverteilung normiert, so ergibt sich 35 



-a'*.[a.in(i-]).(-a*)=M. 



40 



Mit 

45 



50 

folgt: 



55 



Wenn spezieU 5 x = 0 ist, d h, ein Medium mit homogener Leitfahigkeit x vorliegt, vereinfacht sich die 
Potentiaigleichung (7) zur Poisson-Gleichung 

60 

. . 65 

Fur eine verschwindende Strorheinpragung a • J =» 0 wird die Poisison-Gleichung zur Laplace-Gleichung 
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B^<& = 0. (9) 

5 In der praktischen Anwendung der Impedanztomografie wird eine ringformige Anordnung von Elcktroden in 
einer Schnittebene des Untersuchungsobjektes oder ein sogenannter "El ktrodengurter verwendet Die Strom- 
einpragungen und die Spannungsmessungen werden in einer Ebene vorgenonamea Das elektrisch Feld veriauft 
dabei aber nicht nur in der Ebene der Elektroden, sondem erstreckt sich dreidimensional Qber das gesamte 
Volumen des betrachteten Objekts. Die Spannungen an den Elektroden hangen folglich auch von den Leitf ahig- 

10 keitsverhaltnissen im Volumen des Untersuchungsobjektes ab. Eine wesentliche Vereinfachung wird durcb die 
Reduktion der r^umlichen Verhaltnisse auf eine zweidimensionaie Ebene erreicht 

Die Berechnung der elektrischen Feldstarke £ (oder des Potentials O) zu einer gegebenen Verteilung der 
Leitfahigkeit x innerhalb eines Gebietes kann ais "Vorwartsproblem" faezeichnet werdeiL Die Berechnung der 
Verteilung der LeitfEhi^ceit x im Inneren eines Gebietes fur verschiedene nur auf dem Rand bekannte elektri- 

15 sche Fekistarkeverteilungen E (PotentialcUff erenz auf dem Ran<Q stdlt ein sogenanntes Inverses Probtenf dar. 
Das Vorwartsproblem wird durch die direkte Ldsung der Vektorgleichungen {5, 2) 
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65 



9xE=0 



mit Hilfe der Clifford-Algebra der Ebene gelost (vgl. hierzu: Hestenes, D.: New Foundations for Classical 
25 Mechanics (1986X und Hestenes* D, G. Sobczyk: Clifford Algebra to Geometric Calculus (1984)). 
Statt der komplexen Potenttalgleichung (7) 



30 5=0 + (6r)- (a <!>) = — nut R = -ln|— ] 

wird so allgemein die hyperkomplexe Feldgleichung 

35 



(10) 



erhaltea 

Das "inverse Problem" wird zunachst durch eine Linearisierung von Gleichung (10) und eine daran anschlie- 
Bende Interpolation der gemessenen Randwerte von E— Eo in das Innere des Bereichs gelost Danut wird das 
praktische MeBverfahren einer ''Backprojecdon*' entlang der fur homogenes Medium berechneten Aquipotenti- 
45 allinien ahnlich dem Algorithmus von Barber und Brown (vgL oben) — begriindet 

Gegeben sei fur ein ebenes Gebiet mit zwei punktformigen Stromelektroden pi und pz auf dem Rand die 
Gleichung f = Auf punktvektor) 



3,E = (E5^)+^l[5(r-p,)-5(f -pj]. (H) 

Aus kleinen Veranderungen vonR sollen kleine Veranderungen von E — Eo resultieren; 

a,(E-Ej = {Ea,R) = {(E-Eo)a,R>+{Eo a,R>. (12) 



Eine lineare Naherung besteht in der Annahme, daB (E — Eo) 5z R gegenQber (Eo^^, vemachlassigbar ist und 
60 somit naherungsweise 

S,(E-Eo)«(Eoa,R) (13) 



Benotigt wird eine Interpolation Ej der Randwerte yon E — Eo in das Innere der Kreisscheibe, so daB damit 
Gleichung (13) in der Form 
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a,(E-Eo) = (E„a,R) 



(14) 



nach R auf geldst werden kann. 

Der Algorithmus der Backprojection entsteht gerade dann, wenn willkurlich die einfache Interpolation 



E,=MEo(2). M = M(C). C = ln 



Z-P2 



z-p, 



(15) 



benutzt wird, dh. das elektrische Feld Ei entsteht durch eine Multiplikation mh dem Faktor M aus dem 
homogenen Feld Eo. Gleichung(l4) wird dann zu 



10 



15 



Eo(z)aiM=(Eo(z)ajR). (16) 

Nach Anwendung einer konformen Abbildung wird als eine spezielle Losung von Gleichung (16) erhalten: 20 
R=M(C)=f-=^. (17) 

-C-O ^0 25 

Dabei bedeuten E die am Rand des Untersuchungsobjekts gemessene Feldstarke und E© die am Rand des 
Untersuchungsobjektes mit homogener Leitf ahigkeitsverteiiung ermittelte Feldstarke, 

Der Algorithmus der Backprojection erzeugt Bflder von UnterscWeden der Impedanz und nicht Bilder der 
ImpedanzverteDung seibst, d. h. alien Bildem liegt eine Referenzmessung und eine aktuelle Messung zugrunde. 30 
Es gibt zwei Mfiglichkeiten, Unterschiede der Impedanz zu erfassea 

Zum einen konnen Bilder aus zeitlichen Veranderungen, d h.-unterschiedlichen Zustanden gewonnen werden 
(Zustands-differente oder dynamische Abbildungen). Die andere Moglichkeit besteht in der Messung bei rwei 
unterschiedlichen Frequenzen, vorausgesetzt, daB die beiden benutzten Frequenzen zuraindest einen Teil eines 
Dispersionsgebietes einschiieBen (Frequenz-differente Abbildungen). Z.B. zeigt biologisches Gewebe oberhalb 35 
von I kHz strukturbedingt eine ausgepragte Dispersion der Impedanz. Sie ist begrundet in der Existenz zweier 
elektrisch leitender Kompartimente, des extra- und des intrazelhilaren Raimies, die durch elektrisch isolierende, 
sich wie Kapazitaten verhaltende Membranen getrennt sind 

Um einen praktikablen Algorithmus jfur die Biiderzeugung zu entwickein, werden folgende vere'mfachende 
Annahmen gemacht: 40 

( 1) Die Leitf ahigkeitsverteiiung weicht nur wenig von einer homogenen Verteilung ab, 

(2) der Korper hat kreisfdrmigen Umfang und 

(3) die Elektroden sind aquidistant um den Korper plaziert. 

45 

Daruber hinaus konnen bei Stromeinpfagung uber zwei nebeneinanderliegende Elektroden diese bei der 
Berechnung der Aquipotentiallihien zu einem "Stromdipor zusammengef aBt werden, d h. der Abstand zwischen 
den stromfuhrenden Elektroden wird als differentieU klein angenommen. Bei den oben angenommenen 16 
Elektroden sind 16 verschiedene Stromdipole moglicL Diese 16 Oberflacheh-Eiektroden werden in der gc- 
wunschten Schnittebene um den Umfang des Kdrpers aquidistant verteilL 50 

Die Substanz eines Untersuchungsobjektes kann also im einfachsten Fall in bezug auf ihre elektrischen 
Eigenschaften beschrieben werden durch ein Ersatzschaltbild, das aus der Parallelsdiaitung eines Widerstandes 
R und eines Kondensators C besteht Die Impedanz Z betragt dann: 

_ . 55 

1+0 -R-C^ ■ 

Der Quotient von Imaginarteil und Realteil ergibt so 

tanffl = 7-4- = -©RC = Qx = Q — 

Re(2) a . 

mit (p = Phase zwischen eingepragtem Wephselstrom und gemessener Spannung und b = Dielektrizitatskon- 
stante. 
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Wird die Phase <p der Impendanz selbst zur Abbildung herangezogen, erhalt man ein Bild der Verteilung der 
elektrischen Zeitkonstanten RC sowohl bei Zustands-differenten als auch im besonderen bei Frequenz-differen- 
ten Bildern^ wenn bezuglich der Freque nz normiert wird: 



Anstelle der komplexen Impedanz kann auch die komplexe Admittanz ausgewertet werden. Weiterhin ist es 
auch moglich, Kombinationen von RealteiL Imagtnarteil, Betrag oder Phase fur die Abbildung heranzuziehen. 
Die jeweils mit dem Stromgenerator verbundenen Elektroden kdnnen benachbarte Elektroden sein. Es ist aber 
auch mdglich, gegenuberliegende Elektroden fur die Auswertung heranzuziehen. Gleiches gilt auch fiir die 
Messung der Potentialdif ferenzen zwischen den restlichen Elektroden. Es ist mdglich, zwei Elektroden gieidizei- 
tig die Surame von mindestens zwei Strdmen mit unterschiedlicher Frequenz zuzufuhrea Verschiedene Zustan- 
de des Untersuchungsobjektes kdnnen bei einer Frequenz gemessen werden, oder es kann ein Zustand bei 
yerschiedenen Frequenzen gemessen werden. 

Die vom Stromgenerator eingespeisten Strome sind konstant oder durch Messung bekannt Der Stromgene- 
rator kann als Gegentaktstromqaelle ausgebildet sein. Die frequenzabhangigen Eigenschaften des Untersu- 
chungsobjektes sind in einem Querschnittsbild spektroskopisch darstellbar, wobei durch Benutzung von Fre- 
quenzpaaren der gesamte der Messung zugangliche Frequenzbereich zur Spektroskopie ausgenutzt werden 
kann« 

Nachf olgend wird ein Ausf uhrungsbeispiel der Erfindung anhand der Zcichnungen naher erlautert Es zeigen: 
Rg.l eine Impedanztomografie-VorTichtung gemaB einem Ausf uhrungsbeispiel der Erfindung, 
Fig, 2 Aquipotentailinien fur Stromdipole (homogenes Medium), 

Fig. 3 ein mit einem Edelstahlzylinder (0 ^ 20mm) fur eine Frequenz f =40 kHz erhaltenes Bild, 
Fig. 4 ein mit einem Edelstahlzylinder (0 =20mm) und einem Kunststoffzylinder (0=2Omm) fOr f«10 Hz 
erhaltenes Bild, 

Fig. 5 ein mit einem Edelstahlzylinder (0 »^20mm) am Rand erhaltenes Bild, 

Fig. 6 ein Bild eines Thorax, ausgeatmet, bei f =25 kHz, 

Fig. 7 ein Bild ernes Thorax, eingeatmet, bei f = 25 kHz, 

Fig. 8 eine Frequenz-differente Abbildung einer Thorax, 

Fig. 9 eine Frequenz-differente Abbildung eines Unterschenkels, 

Fig. 10 zeitliche Verlaufe eines angewahlten Bildbereiches und 

Fig. 1 1 eine Abwandlung der Auswerteeinheit 

Die Fig. 3 bis 9 zeigen dabei mit der erfindungsgemaBen Vorrichtung gewonnene Bilder. 

Um ein MeB- bzw. Untersuchungsobjekt 17 sind Elektroden 1 bis 16 aquidistant angeordnet 

Wechselstrom wird uber zwei benachbarte Elektroden 1 und 2 in das MeB- bzw. Untersuchungsobjekt 
eingespeist Zwischen den ubrigen 14 Elektroden 3 bis 16 wird jeweils die Differenzspannung zwischen benach- 
barten Elektroden gemessen, also zwischen den Elektroden (3—4), (4—5), (5—6) usw. Es werden also jeweils 
insgesamt 13 Spannungsmessungen ausgefuhrt Dann wird die Stromeinspeisung um eine Elektrode "verscho- 
ben" d h. die Elektroden 2 und 3 sind die Strom zufUhrenden Elektroden. Zwischen den restlichen 14 Elektroden 
werden wieder 13 Differenzspannungen gemessen. Insgesamt wird der Strom nacheinander uber alle 16 mogli- 
chen, nebeneinanderliegenden Elektrodenpaare eingespeist Dadurch wu-d eine Matrix von 16 * 13, dh. 208 
Differenzspannungen ermitteit. 

Bei dieser Vorgehensweise erfolgt die Stromspeisung uber jeweils nebeneinanderliegende Elektroden, d h. 
der Spreizwinkel zwischen den Stromelektroden betragt TUS. Dieser Spreizwinkel ist in Schritten von lUS bis 11 
(diametral angeordnete Stromelektroden) wahlbar. 

Die Vorrichtung enthalt auBerdem eine Ansteuer- und MeBeinheit 18, eine Einheit 19 zur Elektrodenumschal- 
tung, einen Wechselstromgenerator 20, eine Abtast- und Digitalisiereinheit 21 und einen Computer 22 rait 
entsprechendem Programm. Der Computer 22 fuhrt die Steuerung der Vorrichtung, die Berechnung der Impe- 
danz-Komponenten mit Hilfe der Korrelation, die Speicherung der MeBwerte, die Bilderzeugung und die 
Bildverarbeitung durch. Die Bilder selbst werden auif einem Monitor 23 angezeigt 

Der Algorithmus der Ruckprojektion stellt sich bei einem Medium mit inhomogener Leitfahi^keitsverteilung 
in der Praxis folgendermaBen dar (vgL Fig. 2): Zunachst wird der Stromdipol, als Ursache der Aquipotentialli- 
nien, zwischen die Elektroden 1 und 2 gelegt (Fig. 2a). Durch die vom Stromdipol ausgehenden und zu den 
tibrigen Elektroden 3 bis 16 verlaufenden, einem homogenen Meditmi entsprechenden Aquipotentiallinien 
werden "gekrummte Sektoren" definiert Allen Paeln eines solchen Sektors wird diejenige komplexe Leitfahig- 
keit X zugeordnet, die sich aus dem MeBwert der Potentialdif ferenz an den beiden zugehorigen MeB elektroden 
bei homogener Leitfahigkeitsverteilung ergeben wurde, FOr das in Fig. 2a markierte Pixel ist also die zwischen 
den Elektroden 4 und 5 gemessene Potendaldifferenz maBgebend AnschlieBend wird die Stromeinspeisung 
weitergeschaltet; in Fig. 2b zum Beispiel liegt der Stromdipol zwischen den Elektroden 8 und 9. Aus den 
MeBwerten der Potentialdifferenz in alien Qbrigen Elektrodenpaaren ergeben sich wiedenmi x-Werte fur die 
zugehorigen "gekrummten Sektoren^ in einem von ihnen liegt auch das schon in Fig. 2a betrachtete PbceL 

Auf diese Weise wird jedem Pixel P(x,y) fiir jede der n Stromeinspeisungen ein Wert W(npcy) zugeordnet, der 
•entweder der Realteil, der Imaginarteil, der Betrag oder die Phase der komplexen Impedanz sein kann. Unter 
Beriicksichtigung der Gesamtheit aller n Stromeinspeisungen wird dann der Summenwert 



Im(z) 




o Re(z) 
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m ,11 

Q(..y) = ZG(n,x.y). 

als resultierender Pixelwert des Impedanz-Tomogramms definiert Dabei ist 

'W^(n.x.y) -Wg,f (n,x,y) 
WK.f(n,x,y) 
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die relative Differenz der MeBwerte WAkt fur das aktuell untersuchte Medium gegeniiber den McBwerten Wpef 
eines homogenen Referenzmedium und G(ivx,y) ein empirisch festgelegter Wichtungsfaktor mit Werten zwi- 
schen 1 und 10"^, der den Summationsbeitrag der dem Stromeinspeisungsort benachbarten Aufpunkte absenkt is 
Den Werten Q(x,y) wird scMieBlich eine Farb- oder Grauwertskaia zugeordnet, die zwisdien dem Maximal- und 
Minimalwert von Q(x,y) auf gespannt wird. 

Mit der in Fig. 2 gezeigten Vorrichtung wurden zunachst zur Orientierung Messungen an einem elektrolyti- 
schen Trog durchgef uhrt Als Trog dient ein Kunststoffbehalter von 150 mm Durchmesscr and 70 mm Hohe. Um 
den Umfang eines Zyliiiders sind aquidistant 16 Elektroden in einer Hohe von 30 mm fiber dem Boden eingelas- 20 
sen. 

In der Darstellung der tomografischen Bilder kann neben jedes Bild eine Farb-oder Grauwertskaia gesteUt 
werdea An beiden Enden der Skala sind fur die jeweilige MeBgroBe in Prozent die berechneten maximalen 
Unterschiede der aktueilen Messung gegenuber der Referenzmessung angegeben, d. h. die Skala wird fiber den 
ganzen Bereich der berechneten Untersciiiede aufgespannt Die MeBgroBe wird neben den Bildern angegeben 25 
mit \Z\ fQr den Betrag der Impedanz, <p fur die Phase, Re(Z) fur den Realteil und Im(Z) fur den ImaginarteiL 

Die Fig, 3 gibt ein Beispiel f fir ein Zustands-differentes Bild am elektrolytisehen Trog fur einen Edelstahlzylin- 
der (0 =20mm) und einen Kunststoffzylinder (0 «20mm) bei e'mer Frequenz f =40 kHz. Nach der Referenz- 
messung sind der Kunststoffzylinder (A==20nun) in die im Bild rechte Troghilfte und der Edelstahizylinder 
(A « 20mm) in die linke Troghalfte eingebracht worden. Abgebildet sind die Unterschiede in Betrag, Phase, 30 
Realteil und Iraaginarteil der Impedanz gegenuber einer Referenzmessung im nur mit elektrolytischer Losung 
(p = 500 . . lOOOecm) gefiillten Trog. Die MeBfrequenz betrug 40 kHz. Die Konturen der Objekte sind nicht scharf 
wiedergegeben, Um das Zentrum der Objekte erscheinen Zonen jeweils gleichen Widerstandsverhaitnisses, weil 
die Quantisierungsstufen der Farben bzw. Grauwerte verhaltnismaBig groB sind Die hier dargestellten Bilder 
sind geglattet, stets jedoch ungescharft Die angegebenen Prozentzahlen soUten nach der Theorie den Log^ith- 35 
mus der Widerstandsverhaltnisse (aktuelle Messung zu Referenzmessung) wiedergeben. Im verwendeten bilder- 
zeugenden Algorithmus ist aber vorausgesetzt, daB die Referenz eine horaogene Widerstandsverteilung auf- 
wcist und daB die Unterschiede klem bleiben. Folglich ist das Verfahren uberfordert, wenn Widerstandsunter- 
schiede von 6 ... 9 Zehnerpotenzen auftreten. Wenn auch die Prozentangaben hinskhtlich ihrer absoluten 
Betrage mit Vorsicht bewertet werden mussen, so sind sie aber fur .direkte Vergleiche zwischen Biidem sehr 40 
hilfreich. . . • 

Die Fig. 4 zeigt das Ergebnis bei glelcher Anordnung der Objekte im Trog wie in Fig. 3, jedoch mit emer 
MeBfrequenz von 1 0 Hz. 

Bei 40 kHz (Fig. 3) wird der Edelstahlkorper durch den Betrag bzw. den Realteil als Leiter (leitfahiger als der 
umgebende Elektrolyt) korrekt ausgewiesen, durch den Betrag bzw. den Realteil bei 10 Hz jedoch (Fig. 4) als ein 45 
erheblich schlechterer Leiter. Dieser auffallige Befund findet seine Erklarung fiber die Imaginar-Bilder: bei 
beiden MeBfrequenzen sind sie Ihnlich und weisen am Ort des abzubildenden Objekts positive Veranderungen 
gegenuber der elektrolytisehen Losung auf. Die Ursache hierfur liegt in der Grenzschicht zwischen metallischer 
und elektrolytischer Phase begrQndet, die sich wie eine Kapazitat verhalt (kapazitiyer Anteil der Polarisations- 
impedanz). Bei der hohen MeBfrequenz (40 kHz) ist der Scheinwiderstand dieser Kapazitat so niedrig, daB ein 50 
merklicher dielektrischer Verschiebungsstrom durch die Grenzschicht flieBt der sich «n metaliischen Kprper als 
Ladungstragerstrom fortsetzt Deshalb erscheinen beide Komponenten der Impedanz: wegen der metaliischen 
Leitung die reale und wegen der dielektrischen Leitung die Imaginarkomponente. Bei einer niedrigen Frequenz 
(10 Hz) hingegen, wenn der Scheinwiderstand der Grenzschicht hoch ist, wird der Verschiebungsstrom und 
damit der Leitungsstrom durch den Metallkem klein. Folglich werden die Veranderungen des Imaginarteils 55 
klein, und das Realteilbild des Kerns zeigt hoheren Widerstand als der Elektrolyt der Umgebung. 

Im' Phasenbild bei 10 Hz (Fig. 4) erkennt man auf Grund der durch die Kapazitat der Grenzschicht bedingten 
Phasenanderungen nur den Edelstahizylinder. Erst das Imaginarteilbild unterscheidet den Leiter vom Isolator. 
Mit diesen Bildern wird das Phanomen der "Elektrodenpolarisa tion'' veranschaulicht 

Ffir die Gewinnung Frequenz-differenter Bilder ist eine frequenzabhangige Leitfahigkeitsverteilung notwen- 60 
dig. Im elektrolytisehen Trog kann diese Frequenzabhangigkeit durch die Grenzschicht zwischen metallischem 
Leiter und Elektrolyt verursacht werden. Den Biidem von einem Edelstahizylinder (0 -=20mm) am Rand des 
Trogs in Fig. 5 liegenzwei Messungen bei unterschiedlichen Frequenzen (200 Hz und 100 Hz) zu Grunde. 

Bei der Beurteilung der Ergebnisse am elektrolytisehen Trog rauB stets beachtet werden, daB bei Messungen 
in vivo in der Regel weder Edelstahl noch Kunststoffkorper, sondern kleinere Impedanzunterschiede erkannt 65 

werden sollen. ^ u • 

Die Fig. 6 zeigt Bilder aus Betrag, Phase, Real- und Imaginarteii des unteren Ilioraxbereichs (ausgeatmct, bei 
einer Frequenz von 25 kHz). Die Referenzmessung wurde in einem mittleren Atemzustand aufgenommen, d. h. 
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die Lungen sind zum Teil mit Luft gefiillt. Die aktuelle^Messung wird im ausgeatraeten Zustand durchgefOhrt 
Die Bilder stellen den Unt rschied zwischen mittlerer Atemlage und ausgeatmeten Zustand dar, in der die 
Lungenflugel eine positive Lettfahigkeitsanderung aufweisen. Sowohl das Betragsbild als auch im besonderen 
das Phas nbild geben die anatomische Struktur mit zwei Lungenf iugeln gut wieder. 
5 In Fig. 7 ist im direkten Vergleich zu Fig, 6 der eingeatmete Zustand in Bildern von Betrag und Phase ebenf alls 
bei f=25kHz festgehalten. Hier zeigt sich, daB insbesondere die in der Phase enthaltene Information die 
anatomische Struktur gut abbildet, was einen wesentlichen Vorteil der Erfindung darsteDt 

Werden zur Messung der Impedanz eines Objektes zwei verschiedene Frequenzen so gewahlt, daB sie 
wenigstens ein Stuck einer Dispersion einschlieQen, so kann der frequenzbedingte Unterschied der Impedanz 

10 zur Abbildung des Objekts herangezogen werden. Fig. 8 zeigt Bilder aus Frequenz-differenten Messungen am 
Thorax. Diesen Bildem der Schnittebene im unteren Thoraxbereich fiegen jeweils eine Referenzmessung bei 
IQkHz und eine weitere Messung bei 40 kHz zu Grunde. Der Spreizwinkel zwischen den stromfuhrenden 
Elektroden betrug bei der Mefiwertaufnahme dieser Bilder 7-II/a Das Phasenbild gibt besonders gut die 
Struktur der Lungenflugel wieder. 

15 In Fig. 9 rind vom Unterschenkelbereich Frequenz-diff erente Schnittbilder jeweils von Betrag, Phase, Reaheil 
und Imaginarteil der Impedanz dargestellt Das Schienbein liegt im Bereich der Elektroden 13 bis 15 (im 
Uhrzeigersinn, Elektrode 16 ist oben) das Wadenbein in der Nihe yon Elektrode 4. In einem Frequenz-differen- 
ten Bild, ermittelt mit den hier verwendeten Frequenzen von 10 kHz und 20 kHz, werden Knochen nicht 
achtbar, da Kiiochengewebe in diesem Frequenzbereidi keine merkliche Dispersion aufweist. 

20 Besonders in den Bildem von Phase und Imaginarteil werden die anatomisch vorhandenen drei Muskelberei- 
che deudidi getrennt und auf gelost 

Die rait jeder Messung erzeugten Bilder stellen die Unterschiede von zwei Impedanzmessungen in der 
Schnittebene dan In der klinischen Praxis interessiert haufig der zeitliche Verlauf der Impedanzunterschiede 
bestimmter Bereiche des Kdrperquerschnitts. Die Analyse eines ausgewahlten Bildbereichs erfOllt diese Anfor- 

25 derungeningeeigneter Weise. 

Fur eine Detailanalyse kann per Programm ein beliebiger, nicht zwangslaufig zusammenhangender Bildbe- 
reich ausgewahlt werden. In Fig. 10 bt dieser ausgewahlte Bildbereich in den Bildem von Betrag und Phase auf 
der linken Seite weiB markiert 

Im Diagramm wird die Abszisse mit den Nummem der in zeitlicher Folge erzeugten Bilder bezeichnet 

30 Entlang der Ordinate wird die Summe der prozentualen Andemngen jedes Pixels von Betrag, Phase, Realteil 
Oder Emagiiuirteil aufgetragen, dividiert durch die Anzahl der Pixel des gewahlten Bereichs. 

Ein fOr die In^edanztoniografie typischer Anwendungsberekdi ist die Untersuchung der Lungenfunkdon. 
Beim Hnatmen fiillt nichtlettende Luft die LungenfliigeL Entsprechend groB sind dazu die Widerstandsunter- 
schiede im ausgeatmeten Zustand. Die Analyse eines Bildbereichs, in Fig. 10 geeigneterweise uber den Ort 

35 gelegt, an dem die Limgenflugel erwartet werden, macht den zeitlichen Verlauf der Unterschiede von Betrag 
(rot) und Phase (blau) der Impedanz uber eine zeitliche Folge der BDder von zwei Atemzugen deutlich. Die 
Referenzmessung wurde im ausgeatmeten Zustand durchgef uhrt (Zustandsdiff erente Bilder). 

Die Experimente am elektrolytischen Trog und die Messungen in vivo zeigen, daB die Verwendung auch von 
Imaginarteil bzw. Phase der Impedanz, zusammen mit der Moglichkeit, mehrere, relativ niedrige Frequenzen 

40 verwenden zu konnen, zur inhaitlichen Aussage und zur Qualitat der Bilder vorteilhaft beitragen. Es werden 
namlich Strukturen erkannt und getrennt aufgelost, die aus Bildem, die allein aus dem Realteil oder dem Betrag 
der Impedanz erzeugt werden, nicht entnommen werden konnen. 

Trotz aller Schwierigkeiten, die in der Anwendung der Feld- bzw. Potentialtheorie grunden, wird eine 
Auflosung von csl 5% des Durchmessers des betrachteten Querschnitts erreicht Eine Verbesserung der Biidqua- 

45 litat ist durch Weiterentwicklung des bilderzeugenden Algorithmus und durch Anwendung bekannter Methoden 
der Bildverarbdtung zu erwarten. 

Im Vergleich mit den etablierten Tomografieverfahren liegt ein Vorteil der Impedanztomografie in der 
Abbildung mittels bisher nicht genutzter Gewebe-Eigenschaften, namlich der komplexen Leitfahigkeit Weitere 
Vorteile konnen darin gesehen werden, daB das Verfahren relativ preiswert, in der Handhabung einf ach, auch am 

50 Krankenbett anwendbar (leicht transportabel) und — vor allem — nicht invasiv ist 

Fig. 1 1 zeigt ein Ausfuhrungsbeispiel fur die Auswerteeinheit Bei diesem Beispiel werden von den Elektroden 
Ei, E2, ... gewonnene Signale jeweils uber Verstarker, steuerbare Verstarker, Analog/Digital- Wandler und 
Multiplizierer, an denen jeweils ein Referenz-Sinus bzw. Referenz-Cosinus liegt, einem Speicherglied zugefuhrt, 
das der Recheneinheit vorgeschaltet ist Auf diese Weise erhalt die Recheneinheit direkt die Real- und Imaginar- 

55 teile der Impedanz, was fur eine hoheArbeitsgeschwindigkeitsorgt 

In einem digitalen, in der Frequenz von 1 Hz bis 50 kHz programmierbaren Generator werden Sinus- und 
Cosinus-Signale erzeugt Die Sinus-Signale steuem eine Gegentakt-Stromquelle mit Ausgangen I und 1. Die 
Amplitude ist einstellbar von 0 bis 4 mA. Ober Schalter Sa kann der Strom auf zwei beliebige Elektroden 
geschaltet werden, die Polaritat wird durch Schalter Sb festgelegt Die Vorverstarker W dienen als Impedanz- 

60 wandler. Die MeBspannungen aller nicht zur Stromeinspeisung benutzten Elektroden werden uber Schalter Sc 
an die Diff erenzverstarker D V gelegt, in den programmierbaren Verstarkem PV zusatzlich verstarkt und in den 
Analog- Digital- Konvert em ADC digit alisiert Deren Ausgangsdaten werden in den Multiplizierem imd Akku- 
raulierern MA mit den Referenz-Sinus- und Cosinus-Daten korreliert Als Ergebnisse erhalt man jeweils Realteil 
Re(Z) der Impedanz und Imaginarteil Im(Z). Diese werden in dem Speicher abgelegt und zum Computer zur 

65 Bilderzeugung ubertragea 
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* Patenta?lsprucKe 

1. Vorrichtung zum Erzeugen tomografischer Bilder, die die Verteilung der elektrischen Impedanz einer 
Kdrperschicht m einer Ebene (x, y) eines Untersuchungsobjektes (17) darstellen, durch Messung von 
elektrischen Spannungen an der Oberflache des Untersuchungsobjektes (17) bei unterschiedlichen Fre- 5 
quenzen, mit: 

— Elektroden(l — 16), die am Umfangdes Untersuchungsobjektes (17) entlang des Randes der Korper- 
schicht anlegbar sind, 

— einer Ansteuer- und MeBeinheit (18), die Wechselstrom von einem Stromgenerator (20) wenigstens 
zwei der Elektroden (1, 2) zufuhrt, femer zwischen jeweils zwei der restlichen Elektroden infolge des lo 
Wechseistroms vom Stromgenerator auftretende Spannungsctifferenzen miBt und entsprechende MeB- 
signalelieferty 

— einer El^ktrodenumschalteinheit^ die fortiaufend andere Elektroden mit dem Stromgenerator (20) 
verbindet, und 

— einer Auswerteeinheit (22) zum Verarbeiten der MeBsignale ffir eine Anzeige des tomografischen 15 
Bildes auf einer Anzeigeeinheit (23), wobei die Auswerteeinheit (22) jeweiis eine den jeweiligen gemes- 
senen Potentialdifferenzen zugeordnete komplexe Impedanz oder Admittanz mit Real- und Imaginar- 
teil Oder Betrag und Phase ermittelt, 

dadurch gekennzeichnel; 

— daB die Frequenzen so gewihlt sind, daB sie wenigstens ein Stuck des Dispersionsgebiets des 20 
Untersuchungsobjekts einschlieBen, und 

— daB die Auswerteeinheit (22) den Imaginarteil oder statt der Phase das Verh§ltnis von Imaginirteil 
zu Realteil durch die Frequenz dividiert, so daB auf der Anzeigeeinheit (23) ein solches tomografisches 
Bild angezeigt wird, das aus Bildpunkten aufgebaut ist, die durch eine Riickprojektion auf Grund von 
zwei Messungen bei den gewahlten Frequenzen aus mindestens den Imaginirteil und/oder die Phase 25 
enthaitenden Impedanz- oder Admittanzwerten gewonnen sind. 

2. Vorrichtung nach Anspruch 1, dadurch gekennzeichnet, daB der Strom vom Stromgenerator (20) jeweils 
zwei benachbarten Elektroden zugefuhrt wird 

3. Vorrichtung nach Ansprudi 1 oder 2, dadurdi gekennzeichnet» daB die Spannungsdifferenzen zwischen 
untereinander benachbarten Elektroden gemessen werden. 30 

4. Vorrichtung nach einem der Anspriiche 1 bis 3, dadurch gekennzeichnet, daB den Elektroden gieichzeitig 
die Summe von Strdmen mit imterschiedlicher Frequenz zugefuhrt wird 

5. Vorrichtung nach einem der Anspriiche 1 bis 4, dadurch gekennzeichnet, daB jedem Pixel P(x,y) <ies 
anzuzeigenden tomografischen Bildes ein Wert WAkt (n,x,y) zugeordnet ist, der unter BerOcksichtigung der 
Gesaratheit aller N Stromebispeisungen des Summenwertes ^ 35 

^ W^Cn, X, y) - W^^Cn. y) 
Q(x,y) = LGCr.,x,y) w.fc, X, y) 

erhalten ist, wobei 
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W^Cn, X, y) - W^fei, X, y) 45 
W^(n, x, y) ; 

die relative Differenz der Impedanz- oder Adraittanzwerte fur die untersuchte Korperschicht bei zwei der 
gewahlten Frequenzen und G(npc,y) einen empirisch f estgelegten Wichtungsfaktor mit Werten zwischen 1 50 
und 10"^ bedeiiten. 

6. Vorrichtung nach einem der Anspriiche 1 bis 5, dadurch gekennzeichnet, daB die Frequenz des Wechsel- 
stromes zwischen 1 0 Hz und I MHz vorzugsweise unter 50 kHz gewah It wird 

7. Vorrichtung nach einem der Anspdiche I bis 6, dadurch gekennzeichnet, daB die eingespeisten Strome 
konstant oder durch Messung bekannt sind 55 

8. Vorrichtung nach einem der Anspriiche 1 bis 7, dadiux:h gekennzeichnet, daB der Stromgenerator als 
Gegentaktstromquelle ausgebildet ist. 

9. Vorrichtung nach einem der Anspriiche 1 bis 8, dadurch gekennzeichnet, daB die frequenzabhingigen 
Eigenschaften des Untersuchungsobjektes (17) spektroskopisch dargestellt werden. 

10. Vorrichtung nach Anspruch 9, dadurch gekennzeichnet, daB der gesamte der Messung zugangliche 60 
Frequenzbereich zur Spektroskopie ausgenutzt wird 
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